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1^ \ Abstract. We study propagating Alfven waves by solving the time- 

^ . dependent equations of magnetohydrodynamics (MHD) in one dimension 

l/^ \ numerically. In a homogeneous medium the circularly polarized Alfven 

wave is an exact solution of the ideal MHD equations, and therefore it 
does not suffer from any dissipation. A high-amplitude linearly polar- 
ized Alfven wave, on the other hand, steepens and form current sheets, 
\ in which the Poynting flux is lost. In a stratified medium, however, a 

' high-amplitude circularly polarized Alfven wave can also lose a signifi- 

cant fraction of its Poynting flux. 

I ! 
O 

:^ 

o 

c/3 . The acceleration of the fast solar wind is a long-standing problem in solar 

^ ' physics. It is likely that Alfven waves play a significant role in this process 

(e.g. Leer et al. 1982). They can propagate over long distance, which allows 
them to reach the outer corona, because to lowest order they are incompressible 
and do not dissipate. However dissipative damping is required at some point to 
\ avoid too high wind velocities (e.g. Holzer et al. 1983). A closer examination 

shows though that a linearly polarized Alfven wave is compressible to second 
order, since the magnetic pressure, B\/{2^q), varies with half the wave length 
of the Alfven wave itself (Alfven & Falthammar 1963). In a circularly polarized 
Alfven wave, on the other hand, the magnetic pressure is constant along the 
wave, which is the physical reason why the circularly polarized Alfven wave in 
a homogeneous medium is an exact solution to the nonlinear MHD equations. 

Several groups have studied different mechanism through which Alfven 
waves can dissipate. In two dimensions phase mixing due to a transverse gra- 
dient in the phase velocity (e.g. Heyvaerts & Prist 1983) can lead to a strong 
damping of the wave. This mechanism is not available in one dimension though, 
and one has then rather considered the nonlinear coupling of the Alfven wave 
with other modes (e.g. Wentzel (1973), for circularly polarized Alfven waves in 
particular see also Del Zanna et al. (2001) ). Cohen & Kulsrud (1974) found 
in an approximate solution of the nonlinear MHD equations that the linearly 
polarized Alfven wave steepens and can form current sheets. This result has 
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Figure 1. Magnetohydrodynamics waves propagating through a ho- 
mogeneous medium and wave amplitude (^)o = 0.01. In a, b 
Anag = 17 and in c, d /3mag = 0.042. We plot in a and c Ap = p — po 
(solid line) and (dashed line) as a function of distance, z, at 2000 s. 
Ap has been converted to velocity units by multiplying with Cg/po in a 
and with v^/po in c, where po is the background density. In b and d we 
plot Bx (solid line) and Vx (dashed line) as a function of distance z at 
the same time. has been converted to velocity units by multiplying 
with 1/ ^iiqPq. 



later been confirmed by numerical simulations by Boynton &; Torkelsson (1996) 
and Ofman k Davila (1997). 



2. Results 

In this paper wc use the numerical code of Boynton & Torkelsson (1996) (see 
also Torkelsson k, Boynton (1998)) to simulate propagating circularly polarized 
Alfven waves. In the following we will consider both models with a constant 
background density, and stratified models. For all the simulations the temper- 
ature is 10^ K and the period of the waves is 300 s. The isothermal speed of 
sound is 1.2910^ ms . 

2.1. Unstratified models 

For the unstratified models, we differentiate between two cases, when the mag- 
netic field is weak, that is /9mag = 17 and Cg > i'a, and when it is strong that is 
/?mag = 0.042 and vx > Cg. For both the runs p = 1 10~^^ kgm~^. We compare 
the runs in Fig. 1. If we compare our Fig la to Fig. 1 of Boynton & Torkelsson 
(1996) we notice that in both cases there is an acoustic precursor moving at a 
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Figure 2. Evolution of Bx for linearly (a) and circularly (b) polar- 
ized Alfven waves propagating through a homogeneous medium. 



speed close to c^, but in their case the precursor is a sinusoidal oscillation in Ap 
and Vz- Apparently the gradual increase of the magnetic pressure at the front 
of the linearly polarized Alfven wave drives an oscillatory wave, while the sharp 
discontinuity of the magnetic pressure at the front of the circularly polarized 
wave works as a piston. The amplitude of the acoustic precursor in Fig. la 
scales with the square of the amplitude of the Alfven wave. 

When VA > Cs (Fig. Ic and Id) a density enhancement propagates with 
the wave. Vz and Ap are then related hy Vz = va^- While the right edge 
of the enhancement coincides with the front of the Alfven wave, the left edge 
propagates with the velocity Cg. In the leftmost part of Fig. Ic we see the effect 
of the non-linear interaction between the acoustic wave and the Alfven wave. 

Boynton Sz Torkelsson (1996) found that when Bi approaches Bz in strength, 
the Alfven wave starts to steepen into current sheets (Fig. 2a) which enhance the 
dissipation of the wave. The circularly polarized Alfven wave does not develop 
such current sheets (Fig. 2b). 

2.2. Stratified models 

We discus in this paper two models: one with (^)o = 0.1 and one with (^)o = 
1.0, both have a background magnetic field Bz = 3 x 10~^ T. 

At the bottom of the grid va is smaller than Cg and equal to 3.8 x 10'^ ms^^, 
but it increases with z and eventually becomes higher than Cg. Density is taken 
to be equal to 5 x 10~^^ kgm~^ at the bottom of the corona and varies with z 
as 

po(^) = Po(0)exp(-|-^) (1) 

where i^o = 6.1 10^ m is the scale height at z = 0, and R is the solar radius. 

The wave front of the Alfven wave pushes matter upwards forming a density 
perturbation which behaves as an acoustic precursor as long as is smaller 
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Figure 3. Magnctohydrodynamic waves propagating through a strat- 
ified medium with (^)o = 1 and = 3 x 10~^ T. a Ap = p-po (solid 
hnc) and (dashed hnc) as a function of distance, z at time 4375 s. 
Ap has been converted to velocity units by multiplying with Cg/po- b 
Bx (solid line) as a function of distance z at the same time, c The 
z-component of the Poynting flux vector at the same time. 




than Cg. When va becomes larger than Cg the Alfven wave begins to overtake 
the precursor, which however persists and propagates forward with a speed close 
to the speed of sound (Fig. 3a). For a low amplitude Alfven wave the density 
enhancement almost keeps its shape as it propagates while for a high amplitude 
wave the density enhancement changes gradually. 

In Fig. 3c we see a decrease in the Poynting flux in front the acoustic 
precursor. Fig. 4 shows the time average over one Alfven wave period of the 
Poynting flux and the integrated work it does via the Lorentz force. In the 
linearly polarized case the wave loses 70% of its Poynting flux below z = 0.4 R© 
and very little of it is dissipated after that (sec Boynton & Torkelsson (1996) 
Fig. 9) while in the circularly polarized case the wave loses very little of its 
Poynting flux below z = 0.5, but it has lost almost half of it at z = 0.8 R© by 
doing mechanical work on the plasma. 



3. Conclusions 

In constructing models of Alfven wave driven stellar winds it is important to 

understand the mechanism by which the Alfven waves are damped. A linearly 
polarized Alfven wave of high amplitude can steepen and form current sheets 
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Figure 4. Time average over one Alfven wave period of the Poynting 
flux (dashed hne) and the integrated work done via the Lorentz force 
(soHd hne) for a circularly polarized Alfven wave propagating through 
a stratified medium with (^)o = 1 . 

even in a homogeneous medium, which leads to a quick damping of the wave. 
While a circularly polarized Alfven wave does not dissipate its energy in a ho- 
mogeneous medium, we have found that it may do so in a stratified medium. 
The dissipation sets in later than for a linearly polarized Alfven wave though, 
but a significant fraction of its Poynting flux may still be spent on doing work 
on the medium through which the Alfven wave propagates. 
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